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ABSTRACT 

We describe a systematic analysis of energy transfer processes in 
crystals of YAG doped with varying concentrations of Cr and Tm. We use 
both spectral measurements and measurements of the temporal response to 
pulsed excitation to give independent determinations of the microscopic 
interaction parameter for Cr to Tm transfer. We discuss the different 
factors in influencing the temperature dependence of the Cr to Tm 
transfer. Finally, we determine the dependence of the Tm 
cross-relaxation rate on Tm concentration. 

I. INTRODUCTION 

Recently there has been an interest in the development of a solid 
state laser operating in the "eye-safe" region of the spectrum 
(wavelengths greater than 1.54 microns). This interest has led to 
renewed studies of the inter ionic processes among Cr, Tm and Ho ions in 
various crystalline hosts. The ^1^, ^Ig transition of the triply 
ionized Holmium ion is a good candidate for a 2-micron laser as it has a 
large cross-section and a long upper state lifetime. Thulium is expected 
to be a good sensitizer for Holmium because it has an efficient cross 
relaxation and energy transfer coupling to the upper laser level of Ho. 
Furthermore, Thulium absorbs around 780 nm so it could be pumped directly 
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by a semiconductor diode laser or indirectly by energy transfer from 
flashlamp pumped Chromium. In order to obtain pulses short enough for 
several of the proposed applications it is necessary to operate the laser 
in a Q-switched mode. Design of lasers to operate in this mode requires 
a thorough knowledge of the transfer rates for each of the various 
inter ionic processes coupling the energy of the pump to the laser 
transitions. We have initiated a systematic study of the energy transfer 
among Cr, Tm and Ho ions in YAG as a first step in the analysis of 
Q-switched operation of lasers based on the Ho ion sensitized by Cr and 
Tm. In this paper we shall focus on the Cr, Tm system, describing a 
general set of procedures aimed at obtaining detailed information on the 
energy transfer processes. We shall also mention the cross relaxation 
process in Tm. This work provides a clear example of the energy transfer 
processes discussed in references 1 and 2. 

II. CHROMIUM TO THULIUM ENERGY TRANSFER 
II. A. The Transfer Process 

The energy levels of Cr and Tm are shown In Figure 1. We shall be 

interested in the following energy pathway. Energy absorbed in the broad 

4 4 2 2 

Cr absorption bands ( T , T ? ) rapidly goes to the E level. From Cr: E 

1 ^ 3 3 

the energy transfers to the Tm F- and F- levels. From there it 

3 ^ J 

nonradiat ively populates the H^ level. Tm may undergo a cross- 

relaxation ^H^ to the ^F^ level from which it feeds the 

upper laser level in Ho. In this article we shall focus attention on the 
Cr to Tm transfer and the Tm cross-relaxation. 

We are considering a crystalline host in which two Ionic species can 
carry an optical excitation; donor ions D are initially excited and may 
subsequently transfer this excitation to acceptor ions A. The rate of 
energy transfer, by a dipole-dipole transition, from a donor ion to an 
acceptor Ion separated by a distance R is proportional to the Inverse of 
the sixth power of R, [3] 

W_ , = C (6) /R 6 . (1) 

DA 

The coefficient C ^ is called the microscopic interaction parameter. It 
is convenient to express it In terms of the effective decay time of the 
donor Ion (In the absence of any acceptors) and a characteristic 
interaction radius R Q 
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Figure 1. Energy levels of Cr, and Tm in YAG. 
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When the two ions are separated by R Q the rate of energy transfer to an 
acceptor is equal to the decay rate for an isolated donor. 


Generally speaking we are considering two types of measurements - 
steady state spectral measurements and measurements of the temporal 
response to pulsed excitation. We can use each type of measurement to 
get an independent determination of the inter-ionic interaction rate. 


Experiments with pulsed excitation lead to effective decay times of 
the donor ion. This decay time in a sample without acceptors gives the 
lifetime of the donor level and in a sample with acceptors it displays 
the lifetime shortening due to the energy transfer process. Experiments 
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Figure 2. The Tm absorption spectrum and the Cr emission spectrum are 


shown for two temperatures. 





with steady state excitation yield the emission spectra of the donor ions 
and the absorption spectra of the acceptor ion. 


We should mention here that the donor ions may also transfer energy 
among themselves, diffusing the excitation. The more donor diffusion 
competes with donor to acceptor energy transfer the more difficult it is 
to extract information from the examination of the donor luminescence. 
This interference increases with increasing concentration of donors and 
leads to the different regimes discussed in reference 2. 

II. B. Temperature Dependence of Energy Transfer 

An important part of the experimental strategy demonstrated here 
involves examination of the variation of the various physical quantities 
with temperature. There are several ways in which the temperature 
dependence of these phenomena yield valuable information. First of all, 
there is some intrinsic interest in the temperature dependence of these 
factors. In addition, there is evidence that two phenomena may be 
related if they share the same temperature dependence. Finally, a 
statistical enhancement of the data accrues from the examination of the 
tendency of thermal variation. 

II. C. Spectroscopic Measurements 

As discussed in reference 1 there are three separate factors which 
contribute to the microscopic interaction parameter - the oscillator 
strengths for electronic dipole-dipole transitions in the donor ion 
emissions and in acceptor ion absorption and the spectral overlap between 
the donor emission and the acceptor absorption. The first of these, the 
donor oscillator strength, is proportional to the reciprocal of the decay 
time of the donor (in the absence of acceptors); this can be determined 
directly from the decay curves. The second, the acceptor oscillator 
strength Is proportional to the integral of the acceptor absorption 
spectra over the wavelength range of the manifold. The third, the 
overlap integral can be calculated from the donor emission and the 
acceptor absorption spectra. 

Figure 2 shows the Tm absorption spectrum and the Cr emission 
spectrum for two temperatures. While both spectra show broadening with 
temperatures the Cr broadening is much more extensive. There Is a 
sideband developing on the high energy side of the Cr- Emission lines 
which overlaps with the ( 3 H 6 -> 3 H 4 ) Tm absorption. The temperature 
dependence of the overlap integral is due primarily to the development of 
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this sideband in resonance with Tm absorption. The Tm absorption does 
not change significantly with temperature and consequently the Tm 
oscillator strength remains constant with temperature. Using this 
spectral data we can determine the temperature dependence of the Cr-Tm 
overlap and the Tm absorption oscillator strength. The oscillator 
strength of the Cr emission is proportional to the reciprocal lifetime of 
Cr in the absence of Ho. The temperature dependence of each of these 
three factors is shown in Fig. 3. As can be seen there is little change 
of the Tm absorption integral with temperature over the range of 
temperature shown. There is a slight change in the overlap integral at 
lower temperatures reflecting the development of the anti-Stokes 
sideband. At the higher temperatures of interest in laser operation the 
dominant change in the microscopic interaction parameter is due to 
changes in the Cr emission oscillator strength. 

II.D. Time Resolved Spectra 

The time resolved spectra also gives evidence of the energy 
transfer. As can be seen in Figure 4, the temporal response to pulsed 
excitation for Cr alone is exponential. For crystals with higher Tm 
concentration this response becomes non-exponential. At early times it 
displays the rapid transfer of excitation from Cr to those Tm ions which 
are nearby. As time increases and the number of available nearby Tm ions 
decreases, the energy transfer rate diminishes and the decay curve 
settles into an exponential behavior. Evidence for energy migration 
among the Cr ions can be seen in the slopes of the decay curves at longer 
times. As the concentration of Cr increases the magnitude of the slopes 
increases . 

As shown in (1) the time dependence of the decay curve is given by 

ICt) = I p(t) (3) 

where p(t) is given explicitly in reference 1. Taking the natural 

1/2 

logarithm of this we obtain an expression that is quadratic in t and 
depends on three parameters. 

ln(I(t)) = a - (3 \}^ - y t (4) 

The measured data was fit to this expression to find the values of 
(a,/3,y) that minimize mean square error. As shown in reference 1 these 
parameters are 
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Figure 3. The time dependence of the three factors determining the 

microscopic interaction parameter. 1) Cr,Tm overlap integral 
20 5 

x 10 cm , 2) reciprocal lifetime of Cr ( 1/microsec) , 3) Tm 

-3 2 

absorption coefficient integral x 10 (cm ). 



t(psec) 

Figure 4. The temporal response to pulsed excitation for three samples 
of Cr and Tm doped YAG. 






T(K) 


Figure 5. The microscopic interaction parameter determined by both 

spectroscopic measurements and by fitting the decay curves to 
temporal response predicted in Eq.(4). 


a = In ( Iq) 

3 = U/t 0 ) 1/2 (C/C Q ) (5) 

r = + k d 

After obtaining the best fitting, we compare the measured life time with 
the coefficient of t. 

II.E. Cr to Tm Microscopic Interaction Parameter 

As we have discussed, it is possible to determine the microscopic 
interaction parameter by two independent techniques, one using 
spectroscopic data and the other from time resolved date. The 
microscopic interaction parameter calculated from both techniques is 
shown in Figure 5. As evident there is qualitative agreement between the 
results of the two methods. The curved line shown in the figure is a 
best fit to the data points; as discussed previously the temperature 
dependence of the microscopic interaction parameter is due to the three 
independent factors so that the shape of the fitted curve cannot be 
explained further. 


394 






III. THULIUM CROSS-RELAXATION 


Finally, we wish to describe the Tm cross relaxation process. To 
determine the transfer rate we examine the concentration quenching of Tm. 
That is the decrease in lifetime of the Tm level with increasing 
concentration. This concentration quenching effect is shown in Figure 6. 
Fitting the concentration dependence of the lifetime to the quadratic 
f ormula 

T = T 0 / [l + (C/C 0 ) 2 ] (6) 

3 

we can determine the intrinsic lifetime of Tm H„ and the associated half 

4 

concentration (the concentration of Tm at which the intrinsic lifetime is 
decreased by half). Using the intrinsic lifetime we can then separately 
determine the cross-relaxation rate. The concentration dependence of 
this rate Is shown in Figure 7. 


IV. SUMMARY 

In this paper we have described an analysis of energy transfer 
processes occurring in Cr, Tm:YAG. We have used both spectral 
measurements and decay curves to provide separate determinations of the 
microscopic interaction parameter for the Cr to Tm transfer. By 
examining the temperature dependence of the various factors contributing 
to the microscopic interaction parameter we have been able to demonstrate 
that its temperature variation, above 200K, is due primarily to changes 
in the Cr lifetime. Concentration quenching in the Tm cross-relaxation 
led to a determination of the cross-relaxation rate. 

This work was performed at the Solid State Laser Materials 
Laboratory, NASA Langley Research Center, Hampton, VA and supported by 
NASA grants NAG-1-796, NAG-1-955 and NAG-1-957. 
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Abstract 


Recently, research lias been conducted on tlie optical properties of materials as- 
sociated with the development of a solid-state laser in the 2 micron region. In 
support of this effort a mathematical model dcsciibing the energy transfer in a 
holmium laser sensitized with thulium is developed. In this paper we establish 
some qualitative properties of the solution of the model, such as nonnegativity, 
boundedness and integrabilily. A local stability analysis is then performed from 
which conditions for asymptotic stability are obtained. Finally, we repoit on our 
numerical analysis of the system and how it compares with experimental results. 


Introduction 


In the early 1980’s there was a renewal of interest in solid-state lasers, due in 
large part to the development and availability of new host materials. It was around 
this time that NASA began investigating tunable solid-state lasers as promising 
candidates for the Earth Observing System(Eos). During the past several years 
research has been conducted on the sensitized holmium laser, which operates in 
the near infra-red region. It can be used as a source for both LIDAR(L/ght 
Detection And Ranging) and DIAL(D/fTerential Absorption and Lidar) as well as 
on aircraft to make Doppler lidar measurements of windshear since it operates in 
the eye-safe region. 

In a solid-stale laser, a dopant ion substitutes directly into the host lattice. 
When the lanthanide rare earth ion holmium(Ho) is used as a dopant in the host 
crystal yttrium aluminum garnet(Y 3 Als 0 i 2 ), the holmium ion substitutes into 
the yttrium sites. This induces a weak coupling to the host lattice that results 
in narrower absorption and emission features than those commonly observed in 
transition metals. Although the holmium laser has high gain and good energy 
storage properties, ionic interactions among the holmium ions limit the concen- 
tration of holmium possible in any host; therefore to increase the optical energy 
absorbed, a sensitizing ion is included. In the laser system under consideration, 
the lanthanide rare earth ion thulium(Tm) is used as a sensitizer for holmium. 
The interactions between thulium and holmium ions increase the efficiency of the 
Tin-IIo laser but at the cost of introducing more nonlinearities into the model. 


These nonlinearities make both the analysis and the laser dynamics significantly 
different than solid-state laser dynamics studied previously. 

The model comprises thulium ions in energy levels 3 H 4 , 3 F 4 , and 3 H 6 ; holmium 
ions in energy levels 5 I 5 , 5 I 7) and 5 I S ; and the photon density in the optical cavity. 
There are two main types of processes considered in the laser system; namely, 
those processes that are inter-ionic — involving the transfer of excitation from one 
ion to another — and those that are not. Inter-ionic processes include: i. 
relaxation, back transfer, and up-conversion. Processes that are not inter-ionic 
include: spontaneous emission, absorption and stimulated emission. The model 
dbes not account for spatial dependence of excitation in the crystal rod but rather 
utilizes a spacial average over the length of the rod. 

In the next section, the system of equations for the electron populations and 
photon density is introduced and discussed. We then establish some qualitative 
properties of the solution to the system. After that a local stability analysis is 
performed and finally the system is subjected to a numerical treatment and the 
numerical solution compared with experimental results obtained in the laboratory. 

The Model of the Laser Dynamics 


The energy level diagram for the Tm 3+ ion and the Ho 3+ ion in YAG [l] is the 

For 


basis for the idealized model, Figure 1, of lasing action. Insert Figure 1 


the remainder of the discussion the following correspondence is made as a matter 
of notational convenience: thulium energy levels 3 II 6 , 3 F 4) and 3 H 4 correspond 
to energy levels 0,1, and 2, respectively and holmium energy levels 5 I 8 , 5 Iz, and 



5 I 5 correspond to energy levels 0',T, and 2' respectively. The number of thulium 
ions, per cm 3 , in energy level i and at time t will be denoted N,(t) ; z = 0, 1,2. 
Similarly, the number of holmium ions, per cm 3 , in energy level i' and at time t 
will be denoted n,(t) ; ? = 0,1,2 . 

The temporal evolution of the dopant electron populations is described by the 
set of rate equations given by 
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In this system, C, Cj,Ci,7i, and q\ represent the probability of cross relax- 
ation, back transfer, forward transfer, up-conversion, and down-conversion, re- 
spectively. All of these processes are depicted in Figure I. The pumping rate is 
denoted by W p (l) and represents the number of photons, per microsecond, avail- 
able to excite electrons in energy level 0 to energy level 2. 

The quantities 77,77' £ = 1,2 represent the spontaneous emission lifetime of 

energy level 2, A, respectively. The quantity, — , represents the transition rate 

r u 

(due to spontaneous emission) from level i to level j . Similarly, — represents the 

T ij 

spontaneous emission transition rate from level i' to level j ' . These two quantities 


satisfy the relations 


I ^ 1 I (zd 1 n 

- = — and —, = J2— * = 1,2 
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The photon density is denoted by <p{t) and the rate equation describing the 
temporal evolution of the photon density in the optical cavity is given by 


d<t>{t) .. . l ( , .71 , A * 0 , «i(0 

— — = va<f)(t) u,(0 n o (0 + sp 0 . 

(It \ go J T c Tji 


Here, a denotes the transition cross section, v is the velocity of light in the laser 


crystal, <7, is the number of manifolds associated with energy level 1' and likewise 


in 


for go- The contribution clue to spontaneous emission of holmium ions in energy 
level V is denoted sp a and the fluorescent lifetime of the material is denoted tj\. 
The lifetime of a photon in the optical cavity is denoted by r c and can be expressed 
as [2] 

2tJ c 

Tc ~ - In R x R 2 ' 

where i c is the length of the optical cavity, c is the velocity of light in vacuo and 
Ri & R -2 denote the reflectivity of the two mirrors. The electron populations of 
thulium and holmium arc constrained by the relations 

AY = No(l ) + A r i(0 + N 2 {t) and n T = n 0 (t) + ni(t) + n 2 (t) (3) 

’where Nr is the concentration of thulium ions, per cm 3 , and nr is the concentra- 
tion of holmium ions, per cm 3 . 

These “population constraints” can be rewritten as 

Aq(0 = AY - N\(t) - N?(t) and n 0 (i) = n T — nj(f) — n 2 (t) (4) 

Substituting these expressions for A r 0 (i) and n 0 (t) into system (1) reduces the 
system from six to four equations. Then by including the rate equation for the 
photon density (2) we obtain 
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The coupled system (5) is the temporal model of the laser system expressed in 
terms of the physical variables. Wc now normalize system (5) to put it into the 
form used for the analysis and calculations of the subsequent sections. Hence, let 
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be the normalized variables. Rewriting the rate equations (5) in terms of the 
normalized variables (6) yields 
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Qualitative Analysis 


In this section the full set of equations, (7), is modified by excluding two physi- 
cal processes — up-conversion and back transfer. This yields a simplified system for 
which certain qualitative properties of the solution are established. Throughout 
this discussion, the parameters are assumed to satisfy the following (physically 
realistic) relationships : 
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Excluding the above mentioned processes yields the following simplified system of 


rate equations 
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(1-2)+ = 


l-z if \-z>0 
0 otherwise 


( 9 ) 
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The system of equations (8) can be written symbolically as 


y' = f(*,y) 


where 



X 


VV p(0(1 - X - y)+ - — - D x x(l - X - y) + 

12 

y = 

y 

and f(/,y) = 

— — — + 2D[x(l - x - t /)+ - Diy{ 1 - z)+ 
7*21 X] 


z 


+ -D-U/O - -)+ + A(t(1 - z )+ - l)P 

r l 


p 


{&[! - 7(1 - -)+] - ^ j P + ftz 


It now follows directly from [3] that if (i) IV,,^) is continuous and (it) the 
solution vector, y, is bounded then there exists a unique, continuous solution to 
* the system (8). The following four theorems may be viewed as a validation of 
the model since they show that the solution to the system (8) possesses the same 
qualitative properties that one would expect the laser system to possess. The 
proofs of the theorems can be found in [4]. The first theorem establishes the 
nonnegativity of the solution vector y, as well as the ground states of thulium and 
holmium. 

Theorem 1 (Nonnegativity) If \V p (t) > 0 fort > 0 and x(0) > 0, y(0) > 0, 
a:(0) > 0, and P( 0) > 0 then 

(i) If 1 — x(0) — y(0) > 0 then 1 - x(t) — y(£) > 0 for all t > 0; 

(n) If 1 — x(0) — y(0) < 0 then there exists a T such that 1 — x(t) — y(t) > 0 
for all t > T; 
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(iii) (a) If l — x(0) — ?y(0) > 0 then x(t) > 0 and y(t) > 0 for all t > 0 

(b) If 1 — x(0) - y( 0) < 0 then x(t) > 0 and y(t) > 0 for all t > 0; 

(zv) If 1 — z(0) > 0 then 1 — c(/) > 0 for all t > 0; 

(t>) If 1 — z(0) < 0 //zen there exists a T m such that 1 — z(t) > 0 for all t > T* ; 

(tu) z(£) > 0 and P(t) > 0 for all t > 0. 

Theorem 1 stales that the normalized ground levels of thulium and holmium 
remain nonnegative on T = (0,oo) or in other words that the constraints (9) are 
self-enforcing. Consequently, we will drop the “ + ” subscripts in system (8) for 
the remainder of the discussion. 

Theorem 2 (Integrability) Let H /r p (Z) be positive, continuous, and intcgrable 
on I. Then if x( 0) > 0 7 y( 0) > 0, z( 0) > 0, & P( 0) > 0, then x(t), y(t), z(t), & 
P(t) are all intcgrable on T . 

Theorem 3 (Boundedness) If \V p (t) > 0 on 1 then x(t), y(t), z(t), and P(t) 
are all bounded on J . 

Finally, Theorem 4 gives the behaviour of the solution when the pumping term 
W v (t) decays to zero but doesn’t decay quickly enough to be intcgrable on [0,co). 

Theorem 4 Let \V p , x, y , z, and P satisfy the hypotheses of Theorem 1. Fur- 
thermore, let VF p (/) — > 0 as t — * co. Then, the solution vector y — >0 as l — + oo. 
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Stability Analysis 


fc-.i 

L- 


In this section the equilibrium solutions arc obtained and a local stability analysis 
is performed. Before doing this however we make the simplifying assumption that 
the spontaneous emission term, sp Q} has a negligible afFect on the asymptotic be- 
haviour of the system and lienee may be dropped. Furthermore, we will consider 
the case of continuous wave (CW) pumping. These two considerations are tanta- 
mount to taking /? 3 = 0 and IVp(J) = \V p in system (8), where W p is a positive 
constant. Doing this yields the following system of rate equations : 

Y t = ll/ p( 1 - x - v) - ~ - D ' x ( l - 1 _ y) 

-j- = — — +2Dix(l - x - y) - D 2 y(l - z) 

at T')\ r i 

— = 7 + D. x y(\ — z) -f — z) — 1 }P 

dt r, 

f - {Wl-7(l-*)l-^}^ (10) 

By setting the right hand side of equations (10) equal to zero, two equilib- 
rium points are obtained. Equilibrium point 1, denoted I, which lias coordi- 
nates (x 1 ,*/ 1 , z^O) and equilibrium point 2, denoted II, which has coordinates 
, where the coordinates arc expressed in terms of the physical pa- 
rameters of the system and are given in [5]. 
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Having located the equilibrium points we now determine the stability proper- 
ties of I and II. For the stability of I, introduce the new hat variables given by 


m 


m 


1 1 


x = x — i 1 , y = y — y l , z — z — z l , and P = P — 0. Substituting these variables 
into the system of rate equations (10) yields 

^ = ^D 1 x i - W p - - Dj(l - i 1 - y 1 )^ x 

+ (D^'-WJy + D^xf + Dixy 

S ^ / 

9\ 

^ = ^_L + 2j D 1 (l-x- l -y 1 )-2D I x 1 )x 

— — h 2 D\X 1 + £12(1 — c ! y + D2IJ 1 z 

+ D 2 yz — 2D 1 (i) 2 — 2D\iy 

s ^ > 

97 

+ - -jZ X - \)P -D.OJZ - 

> V, ' 

92 

^ = {ft(l -7 + 7^)-^}P + ftT^P 

Pi 

Now, letting .V = [.r, y, £, P) T k g = [yi, y 2 , Oz, Q- i] T , the above system may 
be written more compactly as 

X = AX + g ( 11 ) 

where /l represents the coefTicient matrix of the associated linear system. This is 
„ an “almost linear” system so I will be asymptotically stable if the eigenvalues of A 
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have negative real parts. Invoking the Routh-IIunvitz Theorem [6] the following 
necessary and sufficient condition for I to be asymptotically stable is obtained 

A [ 1 - 7 ( 1 - j ')]--<0 

r c 

The stability analysis for II parallels that of I and yields one condition for II to 
be asymptotically stable. The algebra associated with this condition, however, is 
tremendously complicated and hence we resort to a parametric invest Nation of 
the stability of the equilibrium points. 

The numerics indicate there is an interplay of stability between I and II. To see 
this, consider W p as a parameter and choose it to be “small 11 . By computing the 
equilibrium values of I & II and then checking the Routh-IIurwitz Criteria, the 
stability of each point is determined. By gradually increasing W p and repeating 
the process we see there is an interchange of stability that goes like this. For W p 
small, the P coordinate of II is negative, II is unstable, and I is asymptotically 
stable (i.e. no lasing). As \V p is increased, the P coordinate of II eventually 
* reaches zero and thus the two equilibrium points coalesce (this is the threshold of 
lasing). Finally, .as IF,, is increased still further, the P coordinate of II becomes 
positive, II becomes asymptotically stable, and I becomes unstable (i.e. lasing 
occurs). 

This is illustrated in Figure 2 where P 2 is graphed as a function of W p and the 
stability of I and II is noted. The value of \V p at which the interchange of stability 
occurs is a bifurcation point and is denoted \V * . From numerical calculations it 
was found that W* = 0.342 x 10 -5 . 
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It is also of physical interest to see how the system responds a s the cavity 
lifetime of a photon, r C) is varied. Choosing r c small and following the same 
procedure as above yields the same interplay of stability as before. Figure 3 gives 
P 2 as a function of r c , once again noting the stability of each point. From this we 
see that r c has a bifurcation point which will be denoted r“ and is approximated 


by r c " =4.81 x 10 5 . 


Insert Figures 2 &; 3 


Numerical Analysis 


In this section the numerical solution to system (8) is obtained. This is done 
by employing cither the subroutine LSODA (when using a VAX 11/750) or the 
subroutine DDRIV2 (when using an IBM PC). In either case, all computations 
were performed in double precision arithmetic. LSODA was developed in 1987 at 
Lawrence Livermore National Laboratory in Livermore, California by L.R. Pet- 
zold and A.C. Hindmarsh. DDRIV2 was developed in 1979 and revised in 1987 
by D.K. Kahaner, National Bureau of Standards, and C.D. Sutherland of Los 
Alamos National Laboratory. Both subroutines were created for the numerical 
integration of stiff and nonstilT systems of first order ordinary differential equa- 
tions. When the system is stiff, the subroutines use a backward difference formula 
(BDF) to perform the numerical integration and when the system is not stiff 
they use a higher order Adams method. Both LSODA and DDRIV2 have the 
capability of automatically switching from one method to the other as the sys- 
tem passes from stiff to nonstiff regions. The BDF method was chosen by both 


d.e. solvers throughout the interval of integration which indicates that system 




(S) is stiff. Roughly speaking, a system is stiff when the characteristic equation 
of the associated linear system has a root with a “large” negative real part. A 
more detailed treatment of stiffness can be found in [7]. The numerical values of 
the system parameters are given in Table 1 and were used in the computation of 
the numerical solution. Insert Table 1 Employing the program in [8] the 
numerical solution of (8) was obtained. A typical example of the numerical solu- 
tion is given in Figure 4 where the upper lasing level and the photon u- .. ire 
plotted. Insert Figure 4 Since XV p is constant, a comparison of the “long 
time” behaviour of the numerical solutions with that of the equilibrium solutions 
can be performed. Using the parameter values in Table 1 and the Routh-Hurwitz 
criteria for asymptotic stability it is found that II is asymptotically stable and I 
is unstable. The coordinates of II, as predicted by the stability analysis are given 
in the left column below and the numerical solutions at t = 200/rs are given in 
the right column below. 

.r 2 = 0.1304363 x = 0.1304372 

y 2 = 0.7C3740S y = 0.7637406 

z 2 = 0.523S095 z = 0.5238095 

P 2 = 0.0759492 P = 0.0759491 

From this it is evident that the d.e. solver captures the “long time” behaviour of 
the solutions extremely well. 

Various pumping schemes were considered in the model. Figure 5 gives the 
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numerical solution for P(t) when the pumpin 0 ' term is taken to be 


yV P (t) = a 2 ic~^° 

where a = .01133 . From this figure wc see that the photon density has irregular 
oscillations and decays rather rapidly, becoming negligible by about t — 250/rs. 

Insert Figure 5 

We now consider two modifications that are made to upgrade the mods-! The 
first is to include the terms that account for the back transfer of energy and that 
were excluded at the beginning of the section on qualitative analysis. The second 
is to replace fa & fa in the rate equation for the photon density with 

fc * = ©*■ 

The rationale for the second change is as follows. A mathematical model that 
accurately describes the dynamics of a laser system must account for both spatial 
and temporal variations in the dependent variables. Allowing for both types of 
variation yields a system of nonlinear partial differential equations (p.d.e/). Due 
to the intractability of the equations and the amount of effort and computer time 
involved in solving the system numerically, various averaging techniques are used 
to eliminate the spatial dependence and thus make the system more tractable. L.F. 
Roberts et.al. [9] gives three, such averaging schemes along with their resultant 
o.d.e. models. The conclusion of this work is that only the temporal model 
obtained by taking a spatial average over the optical length of the cavity, C c , yields 
numerical results that agree qualitatively with those predicted by the spatial and 



temporal model. This avcra.!-' 


scheme can ne incorporated into the model by 


H 


t j 

I I 


fc' i 
: * 

u 


I ? 


replacing ft and ft by the quantities given above. Doing this and including the 
back transfer terms yields the system 

~T7 = WpCOO - x - £ ) i.r(l - x - y) 

at T~2 


dv x v 

— = — + 2D t x(l - x - y) - D 2 y{\ - z) + D 3 z(l - x - y) 

at T2 1 T j 

-J 7 = -4 + D. t y{l - z) - D s z(\ - x - y) + ftfrfl - z) - l]P 

at r x 


‘li = {^[1-^0 - ^)1 - + 


( 12 ) 


The numerical solution to the system (12) was obtained giving special attention 

l 

to the qualitative behaviour of the photon density for various values of C[ and — . 
Recall that C' represents the probability of back transfer occurring. The pumping 
term was taken to be of the form 


\V p (t) = a 2 lc~° ‘ . 

In Figure 6 the normalized photon density is graphed for a set of parameter 
values that are “close” to the values used in experimentation. Comparing Figure 6 
with Figure 5 we see a salient disparity in the qualitative behaviour of the photon 
density. The previously observed erratic spiking has been replaced by regular and 
temperate oscillations. Figure 7 is a picture of the energy output of the laser 
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system as displayed on the screen of an oscilK- >cope. 

By comparing Figure 6 with Figure 7 we see that the behaviour of the photon 
density, as predicted by the temporal model (12), is in remarkable agreement with 
the behaviour observed in the laboratory (both having about 7 oscillations in 25 
fis). 

Finally, we further upgrade the model by including the terms associated with 
the up-conversion of energy from the 5 I 7 energy level to the 5 I 5 energy . f 
holmium. This now brings us full circle since we are back to the full set of 
equations (7) with the exception that 0 2 & 0 3 are replaced by 0' 2 &; 0' 3) respectively, 
as-defined above. The numerical solution to the sj'-stem with up-conversion was 
computed taking the probability of up-conversion & down-conversion as 


Insert Figures 6 &: 7 


< 7 , = 5.0 X 1(T 22 & q\ = 1.2 X 1(T 23 


respectively. Figure S shows the interplay between the photon density and the 
upper lasing level of holmium and Figure 9 gives a phase portrait of the same. 


Insert Figures S 


Summary and Conclusions 


We have developed a temporal model of the dynamics of an optically pumped 
co-doped six* level solid state laser. The model was developed to study the inter* 
ionic transfer of energy between thulium and holmium and how this affects the 
performance of the laser system. In this paper we report on the quantitative and 


qualitative behaviour of the solutions. 
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The validity of the simplified model was established by showing that — under 
appropriate conditions on the pumping term — solutions to the system exist and 
are bounded. Further, we showed that if the initial conditions are physically 
realistic and the pumping term is well enough behaved, then the solutions are 
both nonnegative and integrable (integrability is desirable since calculating the 
efficiency of the laser system involves integrating over the photon density). 

Two equilibrium solutions were obtained and a local stability analysis per- 
formed. The interplay of stability between the two equilibrium points was then 
demonstrated parametrically. 

The system was solved numerically and found to be stiff. The d.c. solvers 
used performed exccptio. ily well and gave numerical solutions that agreed with 
the stability analysis and the experimental results. By gradually upgrading the 
simplified model it was found that both back transfer and up-conversion affect 
the time of lasing as well as the magnitude and frequency of the laser output 
but it is an alternate averaging scheme that has the most dramatic afTect on the 
qualitative behaviour of the output pulse. 
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Table 1 System Parameters for Tm-Ho:YAG Laser 




Nt = 1 x 10 21 /cm 3 
n T — 1 x 10 2O /cm 3 
£ = 0.3 cm 

<j>norm = 1 X 10 1 ®/ C? 7 Z 3 

£ c = 17.5cm 
/?! = 1.0 
/?2 = 0.95 


C = 


1 


40 x 10’' Cm3/ " S 
Cl = 475 x 10’° C ’ ,,3/ " S 

c; = o.o 


t-i = 450 /is 
r 20 = 900/is 
r 2i = 900/is 
t i = 11,000/xs 

Tj = 8 , 500 /is 

£o = 1 
9i = 1 
7 = 2 

cr = 7 x 10 -21 cm 2 

7i = 0.0 
9i = 0.0 

v = 30, 000 cm/ /is 
T c = 1 X 10 -3 /iS 
„ = 6 x 10 _3 /cm 3 • fts 

sp 0 = 1 x 10 -G 
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Figure 1 

Idealized model of the Tm-Ho:YAG laser with the following processes represented: 
(a) Stimulated Absorption, (b) Cross Relaxation, (c) Back Transfer, (d) Forward 
Transfer, (e) Up-Conversion, (f) Down-Conversion, and (g) Stimulated Emission. 

Figure 2 

P coordinate of II as a function of the pumping term. 

Figure 3 

P coordinate of II as a function of the cavity lifetime r c . 

Figure 4 

Normalized photon density and upper lasing level with a constant pump. 

Figure 5 

Normalized photon density with a decaying exponential pump. 

Figure 6 

Photon density when using an alternate averaging technique. 

Figure 7 

The energy output of the Tm-Ho:YAG laser as displayed on the screen of an os- 
cilloscope. 

Figure S 

Upper lasing level and photon density with up-conversion included in the model. 
Figure 9 

Phase portrait of the upper lasing level and the photon density. 
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